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Abstract 

Theoretical  investigations  on  the  electronic  structure,  chemical  bonding,  and  ground-state  and  high  pressure  properties  of  metal  and  complex 
hydrides  are  reported.  Recently  a  series  of  metal  hydrides  with  unusually  short  H-H  separations  has  been  synthesized  in  Norway.  The  origin  of 
this  unexpected  short  H-H  separation  has  been  explained  and  several  compounds  with  even  shorter  H-H  separations  have  been  predicted.  From 
a  systematic  study  on  around  100  ZrNiAl-type  compounds  we  have  demonstrated  that  the  electron  localization  function  is  a  powerful  tool  to 
predict  hydrogen  positions  in  metal,  intermetallic  and  alloy  matrices.  Based  on  these  observations  we  have  devised  a  new  empirical  rule,  termed  as 
“site-preference  rule”,  for  prediction  of  suitable  sites  for  hydrogen  accommodation.  First  principle  total-energy  calculations  have  been  carried  out 
for  AH,  EH2,  AEH3,  AXH4,  and  A3XH6  (A  =  alkali  metal,  E  =  alkaline  earth  metal,  X  =  B,  A1  or  Ga)  compounds  to  study  structural  phase  stability 
at  ambient  and  high  pressures.  We  have  reproduced  crystal  structures  of  known  phases  and  predicted  structural  parameters  for  other  compounds 
within  these  series.  We  have  found  that,  if  the  cation  (A+  or  E2+)  radius  is  relatively  small,  one  can  expect  several  pressure-induced  structural 
transitions  with  huge  volume  collapse  at  the  first  phase-transition  point.  These  materials  have  mixed  chemical  bonding  character  which  we  have 
characterized  by  charge  density,  charge  transfer,  electron  localization  function,  site-  and  angular-momentum  projected  density  of  states,  crystal 
orbital  Hamilton  population  analysis,  Born  effective  charges,  and  Mulliken  population  analyses. 
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1.  Introduction 

Hydrogen,  which  has  one  of  the  highest  gravimetric  energy 
densities  of  any  fuel,  unfortunately  is  also  the  lightest  of  all 
elements.  This  means  that  typically  large  volumes  or  high 
pressures  are  required  to  store  appreciable  amount  of  hydrogen 
needed  to  permit  a  fuel  cell  vehicle  to  exhibit  an  acceptable 
operating  range.  Developing  safe  and  reliable  hydrogen  storage 
technology  that  meet  performance  and  cost  requirements  is 
critical  to  use  hydrogen  as  a  fuel  for  both  vehicular  applications 
and  off-board  uses  such  as  for  stationary  power  generation  and 
for  hydrogen  delivery  and  refuelling  infrastructure.  In  order 
to  achieve  [1]  the  storage  capacities  of  2  kWh  kg-1  (6wt.% 
hydrogen)  and  1.5  kWh  L-1  (0.045  kg  hydrogen  L_1)  to  use 
hydrogen  as  a  fuel  for  practical  applications,  the  gravimetric 
and  volumetric  capacities  of  the  hydrogen  storage  materials 
must  clearly  be  increased  than  the  present  levels.  Current 
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on-board  hydrogen  storage  approaches  involve  compressed 
hydrogen  gas  tanks,  liquid  hydrogen  tanks,  metal  hydrides, 
carbon-based  materials/high  surface  area  sorbants  such  as 
metal-organic  frameworks,  and  chemical  hydrogen  storage 
using  complex  hydrides.  Chemical  storage  represents  one  class 
of  potential  solutions  to  the  storage  problem.  Metal  hydrides 
can  release  hydrogen  reversibly  under  acceptable  conditions 
(for  example  LaNis^).  However,  these  materials  are  too  heavy 
and  expensive  for  on-board  vehicle  applications. 

For  the  usage  of  hydrogen  as  an  ideal  fuel  for  many  types 
of  energy  converters  one  requires  efficient  and  safe  storage  of 
hydrogen  in  the  form  of  metal  hydrides.  For  stationary  or  heavy 
duty  vehicle  applications,  traditional  interstitial  metal  hydrides 
have  been  successfully  utilized  [2].  Over  the  past  few  decades,  a 
major  challenge  which  still  remains,  is  to  identify  optimal  inter- 
metallics  for  such  purposes.  Rare-earth  alloys  seem  promising 
owing  to  high  hydrogen  capacity  per  volume  unit  and  their  ability 
to  absorb  hydrogen  under  moderate  conditions  [3,4].  The  hydro¬ 
gen  absorption/desorption  properties  of  these  alloys  are  very 
much  dependent  on  the  constituents  and  the  metal-hydrogen 
bonding  interactions  play  a  major  role  for  the  stability  of  the 
hydrides.  In  order  to  optimize  hydride  phases  for  a  certain 
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application,  an  improved  understanding  of  the  role  of  individual 
alloy  constituents  in  relation  to  hydrogen  and  the  subsequent 
influence  on  electronic  and  structural  properties  are  desirable. 

The  crystal  structures  of  intermetallic  phases  are  often  com¬ 
plex  and  there  are  several  potential  interstices  that  might  accom¬ 
modate  the  hydrogen  depending  on  factors  like  the  size  and 
shape  of  the  interstitial  site,  chemical  nature  of  the  surrounding 
atoms,  and  the  interatomic  distances  to  hydrogen  and  the  coor¬ 
dinating  atoms  [5,6].  Structural  studies  of  intermetallic  hydrides 
have  revealed  that  empirical  rules  can  be  used  to  predict  impor¬ 
tant  features  of  the  hydrogen  sub-lattices  in  a  given  metal  matrix 
[7,8] .  The  searches  for  efficient  hydrogen- storage  metal  hydrides 
have  some  extent  been  hampered  by  the  mental  barriers  which 
empirical  rules  have  put  on  the  thinking.  For  example,  the  radius 
of  the  interstitial  hole  which  hydrogen  is  expected  to  occupy 

o 

should  be  >0.40  A.  A  survey  [8]  of  stable  hydrides  shows  that  the 

o  o 

minimum  H-H  distance  does  not  exceed  2.1  A  (the  “2- A  rule”). 
This  empirical  pattern  is  later  [9]  supported  by  band-structure 
calculations  which  ascribe  the  features  to  repulsive  interaction 
generated  by  the  partially  charged  hydrogen  atoms.  A  practical 
consequence  of  this  repulsive  H-H  interaction  in  metal  hydrides 
is  that  it  puts  a  limit  to  the  amount  of  hydrogen  which  can  be 
accommodated  within  a  given  structural  framework.  So,  if  H-H 

o 

separations  less  than  2  A  would  be  possible  this  could  open  for 
new  prospects  for  potential  candidates  with  higher  hydrogen- 
storage  capacity. 

The  structural  knowledge  for  hydrides  is  by  and  large  very 
limited  owing  to  complexity  in  structural  arrangements  and 
difficulties  involved  in  establishing  hydrogen  positions.  If  one 
know  before  hand  the  knowledge  about  the  site  preference 
of  hydrogen  in  the  host  matrices,  it  will  help  to  design  new 
hydrogen  storage  materials  with  desired  weight  percentage  of 
hydrogen.  The  electron-localization  function  (ELF)  is  a  useful 
tool  to  characterize  chemical  bonding  in  solids.  For  example, 
Savin  et  al.  [10]  have  distinguished  different  kinds  of  bonding  in 
wide  range  of  materials  and  also  demonstrated  that  it  is  possible 
to  identify  the  site  preference  for  hydrogen  in  the  Ca  metal  using 
ELF  analysis.  Using  systematic  ELF  analysis  of  around  hundred 
compounds  in  the  ZrNiAl-type  structural  frameworks  we  have 
confirmed  that  ELF  can  reliably  be  used  to  predict  the  site  pref¬ 
erence  of  hydrogen  in  metals,  intermetallics  and  alloys.  Based 
on  this  observation  we  have  developed  new  semi-empirical  rule 
called  “site-preference  rule”  for  hydrogen  occupation  in  metal 
matrices. 

The  known  reversible  metal  hydrides  which  have  the  desired 
gravimetric  densities,  such  as  MgH2  require  high  temperatures 
to  release  hydrogen.  Currently,  there  is  no  hydrogen  storage 
material  exists  which  can  meet  all  DOE  2010/FreedomCAR 
hydrogen  storage  system  targets.  Recently,  it  has  been  shown 
[11]  that  the  complex  hydride  NaAlH4  can  reversibly  absorb 
hydrogen  at  lower  pressures  and  temperatures  than  MgH2  and 
has  a  higher  gravimetric  capacity  and  lower  cost  than  LaNi5H6. 
Complex  hydrides  form  a  new  class  of  reversible  hydrides  which 
have  not  been  fully  explored.  It  is  possible  that  novel  complex 
hydrides  could  have  improved  properties  for  on-board  hydrogen 
storage.  Complex  hydrides  of  aluminium  have  been  considered 
attractive  as  hydrogen  storage  compounds  due  to  their  large 


hydrogen  content.  Unfortunately,  their  application  in  this  man¬ 
ner  has  been  impractical  as  a  result  of  the  great  difficulties  in 
reversing  the  hydrogen  release  reaction.  Since  workers  in  sev¬ 
eral  laboratories  reported  the  discovery  of  a  number  of  catalysts 
that  improve  the  reversing  of  the  hydrogen  release  in  NaAlH4, 
Na3AlH6,  and  LiAlH4,  interest  for  the  use  of  complex  hydrides 
of  aluminium  as  hydrogen  storage  media  was  rekindled. 

In  the  first  part  of  this  paper  we  present  results  on  metal 
hydrides  which  have  unusually  short  H-H  separation  and  also 
analyzed  the  possible  reasons.  The  predicting  capability  of  den¬ 
sity  functional  theory  on  site  preference  of  hydrogen  in  metal 
matrices  is  also  demonstrated.  The  second  part  deals  with  the 
exploration  of  new  complex  hydrides  from  theoretical  calcu¬ 
lations  and  also  the  structural  stability  of  these  materials  at 
high-pressures  based  on  total  energy  studies.  The  chemical 
bonding  in  selected  complex  hydrides  have  been  analyzed  with 
the  help  of  partial  density  of  states  (DOS),  crystal-orbital  Hamil¬ 
ton  population  (COHP),  charge  density,  charge  difference  and 
electron-localization  function,  Mulliken  population  and  Born 
effective  charges. 

2.  Computational  details 

The  Vienna  ab  initio  simulation  package  (VASP)  [12]  has 
been  used  for  the  total-energy  calculations  to  establish  phase  sta¬ 
bility  and  transition  pressures.  The  generalized  gradient  approx¬ 
imation  (GGA)  [13]  includes  the  effects  of  local  gradients  in  the 
charge  density  for  each  point  in  the  materials  and  which  gener¬ 
ally  gives  better  equilibrium  structural  parameters  than  the  local 
density  approximation  (LDA).  Hence,  we  have  used  GGA  for  all 
our  calculations.  The  ground  state  charge  density  and  energy  are 
calculated  using  a  preconditioned  conjugate  gradient  minimiza¬ 
tion  algorithm  [14]  coupled  with  a  Pulay-like  mixing  scheme 
[15].  Due  to  numerical  instabilities  associated  with  integrating 
the  step-function  character  of  the  0  K  Fermi-Dirac  distribution, 
partial  occupancies  of  the  single-particle  wave  functions  are 
introduced  to  accelerate  the  convergence  [16].  An  energy  cut¬ 
off  of  400  eV  was  used,  in  the  “high  precision”  option  which 
guarantees  that  absolute  energies  are  converged  to  within  a  few 
meV  f.u.-1  To  promote  cancellation  of  errors  when  calculate  the 
total  energy  difference  between  structure,  the  calculations  were 
made  with  same  k  point  density  for  all  structures  within  a  series. 
The  structures  are  fully  relaxed  for  all  volumes  considered  in  the 
present  calculations  using  force  as  well  as  stress  minimization. 
Experimentally  established  structural  data  were  used  as  input 
for  the  calculations  when  available,  or  closely  related  structures 
were  considered  to  identify  ambient  and  high-pressure  phases. 
The  PAW  pseudo-potentials  [17]  were  used  for  all  our  calcula¬ 
tions.  A  criterion  of  at  least  0.01  meV  atom-1  was  placed  on  the 
self-consistent  convergence  of  the  total  energy.  The  present  type 
of  theoretical  approach  has  recently  been  successfully  applied 
[18-21]  to  reproduce  ambient-  and  high-pressure  phases  com¬ 
putationally. 

The  modern  theory  of  macroscopic  polarization  [22]  allows 
the  calculation  of  Born  effective  charge  Z*  using  Berry  phase 
approach  implemented  in  the  VASP  code.  To  have  more  insight 
into  the  bonding  behavior  of  complex  hydrides,  we  have 
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Fig.  1.  Calculated  ELF  in  the  (00  1)  plane  for  LaNiln  and  LaNilnFI  1.333.  Different  crystallographic  site  (4 h  and  3 g)  are  marked  by  arrows.  The  iso-surface  values 
correspond  to  0.6  ELF  value.  Legends  for  different  kinds  of  atoms  are  given  on  the  illustration. 


calculated  Born  effective  charges  using  finite  differences.  The 
chemical  bonding  in  complex  hydrides  are  also  analyzed  with 
the  help  of  Mulliken  population  analyses  as  implemented  in 
CRYSTAL03  package  [23].  For  the  structural  optimization 
for  metal  hydrides  we  have  used  the  full-potential  linear 
muffin-tin-orbital  method  [24] . 

3.  Results  and  discussion 

3.1.  Site  preference  of  hydrogen  in  metal  matrices 

In  order  to  design  new  metal  hydrides  for  hydrogen  stor¬ 
age  applications  a  main  prerequisite  is  the  knowledge  about  the 
site  preference  of  hydrogen  in  the  corresponding  non-hydride 
phases.  The  shortest  H-H  separation  between  two  hydrogen 
atoms  in  a  crystalline  framework  is  limited  by  the  H-to-H  repul¬ 
sion  and  it  is  commonly  believed  that  two  hydrogen  atoms  cannot 

o 

be  located  closer  together  than  some  2.0  A  [8].  This  so-called 

o 

“2- A  rule”  appears  to  hold  for  the  absolute  majority  of  precisely 
described  structures  of  metal  (alloy/intermetallic)  hydrides.  In 
contrast,  a  recent  experimental  study  [25]  on  RNiln  (R  =  La,  Ce, 
Nd;  generally  a  rare-earth  element)  series  shows  that  hydrides 

o 

of  these  compounds  violate  the  “2- A  rule”  vigorously.  Using 
density  functional  calculations  we  have  explained  [26,27]  the 
origin  for  the  violation  and  suggested  that  the  H-H  separation 

o 

may  be  reduced  to  even  below  1.5  A,  by  combinations  of  other 
element.  ELF  is  a  measure  for  the  probability  distribution  of 
paired  electrons,  and  this  function  is  able  to  better  distinguish 
different  bonding  situations  for  electrons  than  the  plain  charge- 
density  distribution  [28].  Savin  et  al.  [10]  demonstrated  that  it 
is  possible  to  identify  the  site  preference  for  hydrogen  in  the  Ca 
metal  using  ELF  analysis. 

Using  the  help  of  ELF,  if  one  able  to  predict  correctly  the  site 
preference  of  hydrogen  in  RNiln  compounds,  then  this  approach 
can  reliably  be  used  to  predict  site  preference  of  hydrogen  in 
metal  matrices.  According  to  experimental  neutron  diffraction 
measurements,  hydrogen  occupies  the  4 h  sites  of  space  group 
P62m  with  full  occupancy  [25].  Interestingly  our  ELF  diagram 


for  LaNiln  shown  in  left  panel  of  Fig.  1  has  higher  value  of  ELF 
at  the  4 h  sites  than  other  interstitial  sites.  The  higher  value  of 
ELF  at  the  4 h  site  indicates  that  electrons  at  that  site  are  hav¬ 
ing  nonbonding  localized  nature  than  other  possible  interstitial 
sites.  It  appears  that  hydrogen  prefers  to  occupy  the  interstitial 
sites  those  have  relatively  more  nonbonding  localized  electrons 
in  metal  matrices.  So,  consistent  with  experimental  observa¬ 
tions  [25],  ELF  analysis  correctly  predicts  the  site  preference  of 
hydrogen  in  LaNiln. 

On  this  background  it  is  interesting  to  visualize  the  ELF  for 
LaNilnHi  333  (see  right  panel  in  Fig.  1)  to  look  for  possible 
interstitial  sites  where  hydrogen  prefers  to  occupy.  With  the  4 h 
site  filled  up  by  hydrogen,  another  interstitial  site  (3g)  with  an 
ELF  value  of  around  0.6  appears  (while  other  interstitial  sites 
still  show  negligible  ELF  levels).  A  recent  study  about  deuter- 
ation  of  LaNiln  under  pressure  (pD2  =  4.6bar)  has  revealed  a 
new  phase,  LaNilnDi  54,  in  which  deuterium  occupies  both  the 
4 h  (96%  filled)  and  3 g  (36%  filled)  sites.  The  structural  opti¬ 
mization  for  the  corresponding  filled-up  LaNiInH2.333  structure 
gave  unit-cell  dimensions  and  positional  parameters  in  satis¬ 
factory  agreement  (viz.  within  the  expected  accuracy  for  DFT 
calculations)  with  the  experimental  findings  for  LaNilnDi  54. 
We  have  also  calculated  the  hydride  formation  energies  for  all 
these  compounds  and  the  conclusions  arrived  from  ELF  analysis 
is  consistent  with  the  results  obtained  from  our  hydride  forma¬ 
tion  energy  calculations  [29]. 

It  was  indeed  the  success  of  ELF  to  correctly  predict  the 
site  preference  of  hydrogen  in  just  discussed  hydride  phases 
led  us  to  undertake  a  more  systematic  study  of  the  ZrNiAl- 
type  family.  Based  on  this  systematic  study  we  proposed  [29] 
a  “site-preference  rule”  for  hydrogen  in  metal  matrices  which 
says  that  hydrogen  prefers  to  occupy  the  interstitial  sites  where 
electrons  have  relatively  more  nonbonding  localized  nature  than 
other  possible  sites  in  metals,  metal  matrices  and  intermetallic 
frameworks.  The  origin  of  the  site-preference  rule  is  that,  when 
hydrogen  atom  is  diffused  into  the  metal  matrices  it  tries  to  find  a 
place  in  the  metal  matrix  where  it  can  find  electrons  that  can  eas¬ 
ily  participate  in  the  bonding  interaction  such  that  it  can  achieve 
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a  more  stable  1  s2  paired  electron  configuration  and  thereby  com¬ 
pletes  its  valence  shell. 

The  insertion  of  hydrogen  in  the  metal  matrix  causes  highly 
anisotropic  lattice  changes;  a  large  expansion  along  c  and  a 
small  contraction  in  the  a  direction.  Among  the  studied  hydrides, 
the  hypothetical  LaPtlnH  1.333  phase  exhibits  the  shortest  H-H 

o 

separation  (1.454  A).  The  optimized  unit-cell  parameters  and 
atomic  coordinates  fit  very  well  with  the  experimental  findings 
for  RNiInHi.333,  R  =  La,  Ce,  and  Nd  [27].  From  the  systematic 
studies  on  series  of  metal  hydrides  we  have  found  that  hydro¬ 
gen  always  accept  charges  from  the  metallic  host  lattice  and 
present  in  the  negatively  charged  condition  at  the  metal  hydrides. 
The  negatively  charged  particle  will  have  repulsive  Coulomb 
interaction  and  hence  the  hydrogen  atoms  in  metal  hydrides  are 
well  separated  from  each  other  (the  H-H  separation  is  gener- 

o 

ally  higher  than  2  A)  compared  to  that  in  hydrogen  molecule 

o 

(the  H-H  separation  is  0.746  A).  In  RNiInHi.333  compounds 
the  electrons  at  the  H  sites  are  polarized  toward  the  electropos¬ 
itive  In  atoms  and  also  the  H-H  repulsive  interaction  has  been 
screened  by  the  rare-earth  ions  present  in  the  systems  [26]. 
These  are  the  reasons  for  the  unusually  short  H-H  separation  in 
RNiInHi.333. 

3.2.  Chemical  bonding  in  complex  hydrides 

3.2.1.  Crystal-orbital  Hamilton-population  analysis 

To  gauge  the  bond  strength  we  have  used  crystal-orbital 
Hamilton-population,  COHP,  analyses,  as  is  implemented  in  the 
TBLMTO-47  package  [35].  The  COHP  (viz.  the  Hamiltonian- 
population-weighted  DOS)  is  identical  with  the  crystal-orbital- 
overlap  population.  Ideally  the  COHP  takes  a  negative  value  for 
bonding  and  positive  for  antibonding  states.  In  order  to  illustrate 
the  bond  strength  between  the  constituents,  the  COHP  between 
atoms  in  KAIH4  is  shown  in  Fig.  2  taking  all  valence  orbitals  into 
account.  From  this  illustration  it  seems  clear  that  the  strongest 
bond  in  KAIH4  is  between  Al  and  H.  There  is  considerable 
covalent  interaction  between  Al  and  H.  Both  bonding  and  anti¬ 
bonding  states  are  present  in  the  COHP  between  the  H  atoms 
at  the  VB,  and  hence  the  bonding  interaction  between  the  H 
atoms  may  be  considered  as  weak.  Owing  to  the  ionic  bonding 
between  K  and  H,  the  COHP  between  these  atoms  is  small.  The 
integrated  COHP  values  give  the  bonding  interactions  between 
atoms  quantitatively  and  those  between  Al-H,  K-H,  H-H,  and 
K-Alare  —0.978,  —0.042,  —0.015,  and  —0.017  eV,  respectively. 
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Fig.  2.  COHP  between  different  constituents  for  KAIH4.  Fermi  level  is  set  to 
zero  and  marked  by  the  vertical  line. 


One  of  the  disadvantages  of  using  complex  hydrides  for  practi¬ 
cal  application  is  the  requirement  of  high  temperatures  to  take 
out  hydrogen  from  these  materials.  From  our  COHP  analysis  it 
is  clear  that  if  one  could  able  to  weaken  the  bonding  interaction 
between  Al  and  hydrogen,  one  can  reduce  the  temperature  to 
take  out  hydrogen. 

3.2.2.  Charge  density,  charge  difference,  and  ELF 

In  order  to  improve  the  understanding  about  the  bonding 
interactions,  the  calculated  valence-charge-density  distribution 
for  KAIH4  is  depicted  in  Fig.  3a.  From  this  illustration  it 
appears  that  Al  and  H  form  an  AIH4  molecule-like  structural 
subunits  and  the  bonding  interaction  between  Al  and  H  is  of 
covalent  character.  Moreover,  Fig.  3a  demonstrates  that  there  is 


Fig.  3.  (a)  Total  charge  density,  (b)  charge  transfer,  and  (c)  electron  localization  function  plot  for  KAIH4. 
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Energy  (eV) 

Fig.  4.  The  calculated  angular-momentum-projected  density  of  states  for 
KAIH4.  The  s-states  are  shown  with  shade  and  the  Fermi  level  is  set  to  zero 
and  marked  by  the  vertical  line. 

no  substantial  charge  density  distributed  between  K  atoms  and 
the  AIH4  structural  subunits.  This  clearly  indicates  the  presence 
of  ionic  bonding.  In  order  to  depict  the  role  of  charge  transfer 
we  have  displayed  the  charge  transfer  plot  (the  difference  in  the 
electron  density  of  the  compound  and  that  of  constituent  atoms 
superimposed  on  the  lattice  grid)  for  KAIH4  in  Fig.  3b.  From 
Fig.  3b  it  is  clear  that  electrons  are  transferred  from  K  and  Al 
to  the  H  sites.  But  the  charge  transfer  from  Al  is  not  isotropic 
as  clearly  seen  from  Fig.  3a.  The  anisotropic  charge  transfer 
from  Al  to  the  H  sites  indicates  the  presence  of  iono-covalent 
bonding  between  Al  and  H. 

The  ELF  is  an  informative  tool  to  distinguish  different  kinds 
of  bonding  interactions  in  solids  [10],  and  the  ELF  for  KAIH4  is 
given  in  Fig.  3c.  The  large  value  of  ELF  at  the  H  site  indicates 
strongly  paired  electrons  with  local  bosonic  character.  The  neg¬ 
ligibly  small  ELF  between  K  and  H,  and  the  small  value  of  ELF 
at  the  K  site  with  spherically  symmetric  distribution  indicate 
that  the  bonding  interaction  between  K  and  H  is  purely  ionic. 
The  ELF  distribution  at  the  H  site  is  not  spherically  symmetric 
and  it  is  polarized  towards  Al  atoms  indicate  the  presence  of 
directional  bonding  between  Al  and  H. 

3.2.3.  Density  of  states  analysis 

The  calculated  partial  density  of  states  for  KAIH4  given  in 
Fig.  4  clearly  shows  that  this  material  is  a  wide  band  gap  insulator 
with  the  band  gap  of  5.2  eV.  Generally  the  band  gap  value  of  all 
the  AAIH4  materials  is  of  the  order  of  5  eV  [19,20].  When  one 
replaces  the  Al  atoms  by  boron  atoms,  the  band  gap  increases 
owing  to  the  increase  of  covalent  bonding  interaction  between 
hydrogen  and  boron.  Hence,  the  band  gap  value  varies  between 
5.5  eV  for  CSBH4  to  7.0  eV  for  LiBH4  [36].  If  one  replaces  the 
Al  atom  with  Ga,  we  found  that  the  band  gap  value  decreases 


Energy  (eV) 

Fig.  5.  The  calculated  total  density  of  states  for  AAIH4.  The  Fermi  level  is  set 
to  zero  and  marked  by  the  vertical  line. 

[37].  The  calculated  band  gap  (Eg)  for  these  compounds  varies 
from  4.57  for  LiGaH4  to  ^5  eV  for  CsGaH4  compound. 

The  valence  band  of  these  compounds  consists  of  two  well 
separated  regions.  In  the  case  of  KAIH4,  the  lower  energy  region 
of  VB  around  —  4  eV  contains  mainly  Al-s  and  H-s  electrons 
with  small  contribution  coming  from  K  atoms.  These  states 
are  the  bonding  states  originating  from  hybridization  interac¬ 
tion  between  Al  and  H.  The  top  of  the  VB  is  dominated  by  H-s 
states  as  well  as  Al-p  states  which  are  energetically  degenerate 
implying  the  presence  of  covalent  bonding  between  Al  and  H. 
Considerably  small  value  of  DOS  at  the  K  site  indicates  the  ionic 
bonding  between  K  and  AIH4. 

The  total  DOS  (Fig.  5)  for  the  AAIH4  phases  comprises  three 
well  separated  regions;  region  I:  below  —5  eV,  region  II:  from 
—3  to  OeV,  and  region  III:  above  5  eV  (unoccupied  states;  see 
Fig.  5).  On  going  from  LiAlH4  to  CSAIH4  the  width  of  the  bands 
(in  particular  the  VB;  see  Fig.  5)  is  narrowed.  Recalling  that  the 
Al-H  distance  does  not  vary  much  along  the  AAIH4  series  it 
seems  safe  to  conclude  that  it  is  the  enhancement  in  the  A-H 
interaction  which  is  the  main  reason  for  the  reduction  in  the  VB 
width.  Also,  on  going  from  Li  to  Cs  systems  the  energy  separa¬ 
tion  between  the  two  peaks  at  the  VB  systematically  increases. 
This  indicates  that  the  ionicity  of  the  system  increases.  The  DOS 
at  the  VB  region  for  LiAlH4  and  NaAlH4  are  having  broad  fea¬ 
tures  than  other  three  compounds  in  Fig.  5  indicating  that  the 
hybridization  interaction  between  atoms  in  these  two  systems  is 
relatively  larger  (Table  1). 
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Table  2 


Mulliken  population  analysis  for  selected  hydrogen-containing  compounds 


Compound 

Atom 

MEC 

Overlap  population 

LiH 

Li 

+0.98 

-0.003  (Li-H) 

H 

-0.98 

MgH2 

Mg 

+1.87 

-0.040  (Mg-H) 

H 

-0.93 

BeH2 

Be 

+1.63 

0.045  (Be-H) 

H 

-0.82 

A1H3 

Al 

+2.22 

0.091  (Al-H) 

H 

-0.74 

LiAlH4 

Li 

+1.01 

0.171  (Al-H) 

Al 

+2.01 

-0.021  (Li-H) 

H 

-0.75 

Li3AlH6 

Li 

+1.01 

0.105  (Al-H) 

Al 

+2.08 

-0.020  (Li-H) 

H 

-0.85 

The  Mulliken  effective  charges  (MEC)  are  given  in  terms  of  e. 


3.2.4.  Mulliken  population  analysis  and  Born  effective 
charges 

In  order  to  make  a  quantitative  conclusion  about  bonding 
interactions  between  constituents  it  would  be  useful  to  identify 
the  amount  of  electrons  on  a  particular  atom  and  populations 
between  atoms.  Although  there  is  no  unique  definition  of  how 
many  electrons  are  associated  with  an  atom  in  a  molecule  or  a 
sub-unit  of  a  solid,  it  has  nevertheless  proved  useful  in  many 
cases  to  perform  population  analyses.  Due  to  its  simplicity  the 
Mulliken  [38]  population  scheme  has  become  the  most  familiar 
approach  to  count  electrons  associated  with  a  given  atom.  How¬ 
ever,  this  method  is  more  qualitative  than  quantitative,  giving 
results  that  are  sensitive  to  the  atomic  basis.  Mulliken  charges 
are  reported  in  Table  2  for  some  saline-like  (LiH,  BeH2  and 
MgH2)  and  complex  hydrides.  LiH  is  considered  as  a  pure  ionic 
compound  and  the  calculated  Mulliken  charges  reflect  nearly 
pure  ionic  picture  with  Li+  and  H-.  The  overlap  population 
between  Li+  and  H-  is  also  close  to  zero,  as  expected  for  an 
ionic  compound.  Similarly,  the  case  of  MgH2,  BeH2,  and  AIH3 
the  bonding  interaction  is  mainly  ionic,  but  the  degree  of  ionicity 
is  reduced  from  MgH2  to  BeH2  and  further  to  AIH3,  viz.  these 
compounds  exhibit  some  covalent  character  as  evidenced  by  the 
non-zero  overlap  population.  The  Mulliken  effective  charges  for 
Li,  Al,  and  H  in  LiAlH4  and  Li3  AlH^  indicate  that  the  interac¬ 
tion  between  the  Li  and  AIH4/AIH6  is  ionic.  There  is  a  finite 
overlap  population  between  Al  and  H  within  the  AIH4/AIH6 
units  which  reflects  the  presence  of  partial  covalent  character 
of  the  Al-H  bond.  However,  the  magnitude  of  overlap  popula¬ 
tion  is  smaller  than  that  for  pure  covalent  compounds.  Also,  the 
partial  charges  (around  two  electrons  transferred  from  Al  to  H) 
imply  the  presence  of  significant  ionic  contribution  in  the  Al-H 
bond. 

Born  effective  charge  is  another  quantity  to  analyze  the 
bonding  character  of  the  materials.  We  used  the  King-Smith 
and  Vanderbilt  [39]  method  to  calculate  the  polarizations  of  per¬ 
turbed  cells  and  from  this  the  Born  effective  charge  tensors  of  the 
ions  have  been  calculated.  In  Table  3  the  calculated  Born  effec- 
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Table  3 


Calculated  Born  effective  charge  tensors  for  Li3  AITC 


yy 

zz 

xy 

yz 

zx 

xz 

zy 

zLi* 

0.990 

0.957 

0.985 

0.017 

0.102 

0.014 

0.061 

0.024 

0.102 

Zaii* 

2.237 

2.237 

1.735 

0.034 

-0.001 

0.000 

-0.001 

0.000 

-0.001 

Za12* 

2.024 

2.024 

2.251 

0.024 

-0.001 

0.000 

0.000 

0.000 

-0.001 

* 

X 

N 

-1.018 

-0.738 

-0.846 

0.082 

0.027 

-0.174 

-0.190 

0.073 

0.027 

ZH2* 

-0.911 

-0.831 

-0.817 

-0.088 

0.125 

0.210 

0.153 

0.116 

0.125 

five  charge  tensors  for  the  constituents  in  Li3  AlH^  are  tabulated. 
For  Li  and  Al,  the  diagonal  components  of  effective  charges 
in  the  Cartesian  frame  are  Zxx  =  Zyy  -=f  Zzz  and  the  off-diagonal 
components  are  negligible.  However,  at  the  H  site,  finite  value 
of  effective  charges  are  present  in  the  diagonal  as  well  as  off- 
diagonal  components  along  with  their  anisotropic  distribution 
implies  that  a  finite  directional  bonding  is  present  in  this  mate¬ 
rial.  If  the  interaction  is  purely  ionic,  then  one  can  expect  that  the 
effective  charges  along  the  diagonal  component  should  be  equal 
and  the  off-diagonal  may  vanish.  It  can  be  seen  from  Table  3  that 
the  dynamical  charge  for  Li  is  close  to  the  nominal  ionic  charge 
(+1  for  Li).  For  Al,  its  dynamical  charges  are  1.735-2.24, 
which  are  smaller  than  the  nominal  charge  of  +3.  This  indicates 
that  part  of  the  Al  charges  also  participate  in  the  directional 
bonding. 

3.3.  Prediction  of  new  complex  hydrides 

Considerable  research  is  being  conducted  to  develop  hydride 
compounds  that  can  store  hydrogen  with  high  volumetric  and 
gravimetric  densities,  and  charge/discharge  it  reversibly  under 
the  temperature  of  80-100  °C  and  1  bar  absolute  pressure  oper¬ 
ation  conditions.  The  current  renaissance  in  hydrogen  storage 
research  was  sparked  by  the  discovery  of  catalyzed  NaAlH4 
compositions  that  have  theoretically  5.6  wt.%  accessible  hydro¬ 
gen  and  are  reversible  under  moderate  temperature  and  pres¬ 
sure  conditions  [12].  Recent  experimental  evidence  shows  that 
reversible  hydrogen  absorption/desorption  proceeds  smoothly  in 
KAIH4  without  introduction  of  a  catalyst  [40].  From  this  obser¬ 
vation  it  appears  that  it  may  be  possible  to  have  complex  hydrides 
with  good  kinetics  those  are  yet  to  be  identified.  In  this  back¬ 
ground  we  are  exploring  new  complex  hydrides  using  density 
functional  calculations.  So  far  we  have  searched  new  compounds 
in  AXH4  (A  =  Li,  Na,  K,  Rb  or  Cs;  X  =  B,  Al,  or  Ga),  A3XH6, 
and  EXH5  (E  =  Be,  Mg,  Ca,  Sr  or  Ba)  series.  From  this  theoreti¬ 
cal  exploration  we  have  identified  several  new  hydrogen  storage 
materials  and  their  ground  state  as  well  as  high-pressure  struc¬ 
tural  properties  are  predicted. 

The  uncertainty  associated  with  the  determination  of  co¬ 
ordinates  of  hydrogen  by  X-ray  diffraction  is  large,  and  the  thus 
derived  metal-hydrogen  distances  tend  to  come  out  too  short 
or  deviate  considerably  from  the  reality.  In  general,  compared 
with  systems  such  as  oxides,  structural  understanding  of 
hydrides  is  very  defective  and  fragmentary  [41].  It  is  now  well 
documented  that  first  principle  calculations  can  reliably  be  used 
to  predict  structural  parameters  of  hydrogen  storage  materials 


[19-22,27,36,37,42-44].  From  the  systematic  studies  on  struc¬ 
tural  phase  stability  of  around  50  compounds  we  found  that  these 
compounds  generally  have  tetrahedra  (XH4)  or  octahedra  (XH6) 
structural  subunits.  Depending  upon  the  radius  of  the  cations 
these  polyhedra  get  tilted  and/or  rotated  which  make  different 
crystal  structures  for  these  compounds.  The  calculated  total 
energy  versus  volume  curves  for  RbAlH4  and  CSAIH4  are  given 
in  Fig.  6a  and  b,  respectively.  The  structural  parameters  for  these 
two  compounds  are  not  available  for  comparison.  Interestingly 
both  these  compounds  stabilize  in  the  orthorhombic  structure. 
The  structural  parameters  for  all  the  AAIH4  compounds 
obtained  from  our  theoretical  explorations  are  given  in  Table  4. 
For  LiAlH4,  NaAlH4,  and  KAIH4  the  structural  parameters  are 
available  experimentally.  From  the  table  it  should  be  noted  that 
the  deviation  in  the  predicted  structural  parameters  from  exper¬ 
imental  values  for  these  compounds  are  within  3%.  In  the  case 
of  KAIH4  we  have  first  predicted  [43]  the  structural  parameters 
and  followed  by  our  report,  neutron  diffraction  study  [34]  has 
been  made  on  this  compound.  The  structural  parameters  for 
RbAlH4  and  CSAIH4  are  given  in  Table  4  and  no  experimental 
results  are  available  for  these  compounds.  We  hope  that  the 
present  results  will  motivate  experimental  studies  on  these 
systems. 

Both  RbAlH4  and  CSAIH4  crystallize  in  the  orthorhombic 
KGaH4-type  structure  with  space  group  Pnma.  This  structure 
consists  of  slightly  distorted  [A1H4]-  tetrahedra  which  are  sep¬ 
arated  by  intermittent  Rb/Cs  cations.  In  KAIH4,  the  interatomic 
Al-H  distance  within  the  [AIH4]-  tetrahedra  varies  only  by 
0.9%  and  the  bond  length  is  as  expected  for  an  anionic  com¬ 
plex  of  Al  and  H  [45].  Each  K+  cation  is  surrounded  by  12  H 

o 

atoms  at  distances  varying  between  2.717  and  3.204  A.  In  the 
case  of  RbAlH4,  we  found  that  the  energy  difference  between 
the  ground  state  structure  and  a-LiAlH4-type  structure  (Fig.  6) 
is  very  small.  We  have  identified  a  high-pressure  modification 
around  7  GPa  in  a-RbAlH4  and  it  transforms  in  to  (3-RbAlH4 
(|3-NaAlH4-type  structure)  and  the  volume  collapse  at  the  phase 
transition  point  is  relatively  larger  value  of  12.6%  (see  Table  4 
and  also  Fig.  7). 

It  is  interesting  to  note  that:  at  2.6  GPa  a-LiAlH4  (proto¬ 
type  structure)  transforms  to  (3-LiAlH4  (a-NaAlH4-type)  and 
a  subsequent  transition  from  (3-  to  y-LiAlH4  (KGaH4-type)  is 
established  at  33.8  GPa.  The  important  aspect  of  the  a-to-(3  tran¬ 
sition  in  LiAlH4  is  the  small  energy  difference  [19]  between 
the  involved  phases  (only  11.56  me V  f.u.-1),  or  1.154kJmol_1, 
which  is  indeed  much  smaller  than  that  found  for  other  hydrides 
[21]  and  oxides  [46].  At  the  a-to-(3  transition  point  for  LiAfirG 
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Table  4 

Calculated  cell  parameters  (A),  unit-cell  volume  (A3),  positional  parameters,  bulk  modulus  (Bo  in  GPa)  and  pressure  derivative  (#q),  and  volume  discontinuities  (%) 
at  phase  transition  points  for  MAIH4  compounds 

Compounds  (structure  Unit-cell  dimension  Unit-cell  volume  Positional  parameters  Bo  B'{)  AV 

type;  space  group) 


a-LiAlH4  (proto  type; 
P2\/c) 


a  =  4.8535  (4.8174)a  65.24 

b  =  7.8259  (7.8020)a 
c  =  7.8419  (7.8214)a 
>8=  111.878°  (112.228°)a 


P-UAIH4  <i  =  4.6611b  (4.7312)c 

(a-NaAlH4-type; 

I4\/a) 

c=  10.5219b  (10.7161)c 


y-LiAlH4  (KGaH4-type; 
Pnma ) 


a  =  5.4421 


b  =  4.4843 


59.96 


134.77 


c  =  5.5225 


a-NaAlH4  (proto  type; 
I4\/a) 


a= 4.9965  (4.980  l)d 
c  =  11.0828  (11 . 1483)d 


69.15 


P-NaAlH4 

(SrMgH4-type;  Cmc22 ) 


a  =  3.5493 
13.8304 


62.75 


c  =  5.1133 


KA1H4  (KGaH4-type; 
Pnma ) 


a  =  8.8249  (8.814/ 


/?  =  5.8590  (5.8 19)e 


95.49 


c  =  7.3872  (7.33 l)e 


a-RbAlH4  (KGaH4-type;  a  =  9.4708  107.39 

Pnma ) 

b  =  5.9075 
c  =  7.6779 


(3-RbAlH4  a  =  3.9933  94.35 

(SrMgH4-type;  Cmc22 ) 


0:0.5699,0.4652,0.8245  (0.5603,  12.95  4.10 

0.4656,  0.8266)a 

Al:  0.1381,  0.2017,  0.9319  (0.1386, 

0.2033,  0.9302)a 

HI:  0.1807,  0.0986,  0.7630  (0.1826, 

0.0958,  0.7630)a 

H2:  0.3542,  0.3723,  0.9777  (0.3524, 

0.3713,  0.9749)a 

H3:  0.2361,  0.0810,  0.1146  (0.2425, 

0.0806,  0.1 148)a 

H4:  0.7948,  0.2633,  0.8717  (0.7994, 

0.2649,  0.8724)a 

Li:  (0,  1/4,  5/8)b,c;  Al:  (0,  1/4,  l/8)b’c  25.64  4.35  17 


H:  (0.2527,  0.4237,  0.5413)b  ( 

0.2492,  0.4191,  0.5429)c 

Li:  (0.2428,  1/4, 0.2467)b  (0.2497,  14.25  4.85  4 

1/4,  0.2502)c 

Al:  (0.5120,  1/4,  0.8221)b  (0.5002, 

1/4,  0.7361)c 

HI:  (0.3067,  1/4,  0.9617)b  (0.2815, 

1/4,  0.9617)c 

H2:  (0.7162,  1/4,  0.9631)b  (0.7189, 

1/4,  0.9467)c 

H3:  (0.4889,  0.9833,  0.2943)b 
(0.4998,  0.9173,  0.3279)c 

Na(4a):  0,  1/4,  1/8  19.29  3.9 

Al(4b):  0,  1/4,  5/8 

H(16f):  0.2199,  0.3710,  0.5639 

(0.2372,  0.3869,  0.5456)d 

Na(4a):  0,0.1695, 0.2409  38.41  3.1  4 


Al(4a):  0,  0.4139,  0.2047;  Hl(4a):  0, 

0.3476,  0.4877 

H2(4a):  0,  0.3105,  0.0116;  H3(4a):  0, 

0.0751,0.6811 
H4(4a):  0,  0.5269,  0.3739 

K(4c):  (0.1778,  1/4, 0.1621)b  10.34  4.61 

(0.1775,  1/4,  0.1598)f 

Al(4c):  (0.5663,  1/4,  0.8184)b 

(0.5659,  1/4,  0.8201)f 

HI (4c):  (0.4034,  1/4,  0.9184)b 

(0.4063,  1/4,  0.961  l)f 

H2(4c):  (0.7055,  1/4,  0.9623)b 

(0.7153,  1/4,  0.961  l)f 

H3(8d):  (0.4194,  0.9810,  0.3127)b 

(0.4181,0.9701,0.3137/ 

Rb(4c):  0.1816,  1/4,0.1597  9.2  5.2 

Al(4c):  0.0615,  1/4,  0.6863 
HI (4c):  0.0983,  3/4,  0.3989 
H2(4c):  0.3117,  3/4,  0.0389 
H3(8d):  0.0802,  0.0238,  0.8121 

Rb(4a):  1/2, 0.6206, 0.2833  14.1  3.5  12.6 


b=  14.6472 


Al(4a):  1/2,  0.1154,  0.7607 
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Table  4  ( Continued ) 


Compounds  (structure 
type;  space  group) 

Unit-cell  dimension 

Unit-cell  volume 

Positional  parameters 

Bo 

Bo 

AV 

c  =  6.4933 

HI (4a):  1/2,  0.7996,  0.0670 

H2(4a):  1/2,  0.1717,  0.9990 

H3(4a):  1/2,  0.5992,  0.7814 

H4(4a):  1/2,  0.9888,0.1074 

a-CsAlH4  (KGaH4-type; 

a  =  9.9249 

123.4 

Cs(4c):  0.1847,  1/4,  0.1629 

8.6 

4.6 

3.8 

Pnma ) 

£>  =  6.1988 

Al(4c):  0.0594,  1/4,  0.6970 

c  =  8.0198 

HI (4c):  0.1087,  3/4,  0.3821 

H2(4c):  0.3125,  3/4,  0.0418 

H3(8d):  0.0812,  0.0180,  0.8256 

(3-CsA1H4 

<3  =  4.3887 

108.1 

Cs(4a):  0,  3/4,  1/8 

9.5 

5.1 

(a-NaAlH4-type) 

b=  11.2137 

Al(4b):  1/2,  3/4,  7/8;  H(16h):  0.7622, 

3/4,  0.9434 

y-CsAlH4  (SrMgH4  type; 

<3  =  3.6619 

85.3 

Cs(4a):  1/2,  0.6159,  0.2871 

9.1 

4.9 

9.3 

Cmc22 ) 

b=  14.2468 

Al(4a):  1/2,  0.1221,0.7546 

c  =  6.5 1766 

HI (4a):  1/2,  0.79419,  0.08073 

H2(4a):  1/2,  0.1649,  0.9874 

H3(4a):  1/2,  0.6072,  0.7706 

H4(4a):  1/2,  0.9789,  0.1186 

a  Experimental  values  from  Ref.  [31]. 
b  Calculated  values  at  equilibrium. 
c  Calculated  values  at  transition  point. 
d  Experimental  values  from  Ref.  [32]. 
e  From  XRD  measurements  [33]. 
f  From  neutron  diffraction  [34] . 


the  estimated  difference  in  cell  volume  is  ca.  17%  (Fig.  7)  [19]. 
The  Raman  scattering  measurement  on  the  high-pressure  phase 
also  indicates  the  presence  of  a  large  volume  collapse  at  the 
transition  point  [47].  Comparison  on  the  basis  of  the  equilib¬ 
rium  volumes  for  a-  and  (3-LiAlH4  show  an  even  larger  volume 
difference  (viz.  a  huge  value  of  22%),  which  to  our  knowledge 
is  the  first  establishment  of  such  huge  volume  change  at  a  tran¬ 
sition  point  in  hydrides.  For  example,  in  NaAlFLt  the  calculated 
volume  difference  at  the  transition  point  between  a  and  (3  phases 
is  less  than  4%  [20].  It  should  be  noted  that  (3-  and  y-LiAlFLt 
have  almost  the  same  volume  at  the  (3-to-y  transition  point.  The 
relatively  small  equilibrium  volume  of  p-LiAlFC  along  with  its 
high  weight  content  of  hydrogen  implies  an  increased  hydrogen 
storage  capacity  and  therefore  it  would  be  of  interest  to  explore 
the  possibility  of  stabilizing  this  phase  by  chemical  means. 
Moreover,  this  phase  may  have  improved  kinetics,  because  the 
bonding  behavior  of  the  (3  phase  is  drastically  different  from  that 
of  the  a  phase. 

We  have  calculated  the  total  energy  as  a  function  of  unit¬ 
cell  volume  for  12  different  volumes  and  fitted  the  results  to 
the  so-called  universal  equation  state  [48]  and  thus  obtained 
bulk  moduli  (Bq)  and  its  pressure  derivatives  (By)  for  all  the 
AAIH4  compounds  at  ambient  and  high  pressures  (see  Table  4). 
All  these  compounds  exhibit  very  small  bulk  modulus  at  the 
ground  state  structure  compared  with  metal  hydrides  implying 
that  AAIH4  are  very  soft  and  easily  compressible  materials. 
The  co-ordination  number  (CN)  of  Al  changes  from  4  in  a- 


RbAlH4  to  6  in  the  (3-phase,  viz.  consistent  with  the  expected 
general  trend  of  increasing  CN  with  weakening  of  the  cova¬ 
lent  bonding  by  increasing  pressure.  Bq  for  p-RbAlFC  is  almost 
twice  that  of  the  a-phase  and  this  is  associated  with  large  vol¬ 
ume  collapse  at  the  phase  transition  point.  Also,  the  slope  of 
the  pressure-volume  curve  (see  Fig.  7)  for  (3-RbAlFLt  is  larger 
than  that  of  a-RbAlFLt  indicating  that  the  high-pressure  phase  is 
intrinsically  harder  than  the  ground  state  phase  and  this  is  largely 
due  to  the  increased  CN  and  the  lower  equilibrium  volume 
of  (3-RbAlH4. 

Like  a-RbAlFU,  a-CsAlFU  also  stabilizes  in  the  orthorhom¬ 
bic  structure  with  space  group  Pnma.  The  bulk  modulus  for 
a-Cs AIH4  is  smaller  than  all  the  other  AAIH4  compounds  which 
is  associated  with  the  large  equilibrium  volume.  Moreover,  the 
ionic  radius  for  Cs  is  larger  than  other  alkali  metals  and  also 
the  valence  electrons  of  Cs  are  more  delocalized  than  other 
alkali  metals.  These  make  a-CsAlFLt  the  softest  among  all  the 
AAIH4  compounds.  Similar  to  LiAlFLt,  the  total  energy  stud¬ 
ies  (see  Fig.  6b)  on  CsAIFU  also  show  two  structural  phase 
transitions  one  at  4  GPa  from  orthorhombic  a-CsAlFLt  to  tetrag¬ 
onal  |3-CsA1H4,  and  around  8  GPa  from  tetragonal  |3-CsA1H4 
to  orthorhombic  y-CsAlFU  with  the  volume  collapse  of  3.8 
and  9.3%,  respectively.  Unlike  other  AAIH4  compounds  the 
bulk  modulus  for  CsAlFLt  in  all  the  three  phases  are  almost 
the  same. 

Similar  to  alkali  aluminohydrides,  alkaline-earth  hydrides 
with  low  cation  radii  also  exhibit  a  huge  volume  collapse  at 
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Fig.  7.  Equation  of  states  data  for  alkali  aluminohydrides. 

the  first  pressure-induced  phase  transition  point.  Application  of 
pressure  transforms  a-  to  |3-BeH2  at  7.07  GPa  with  a  calculated 
volume  discontinuity  at  the  transition  point  of  ca.  19.2%  (ca. 
26.4%  equilibrium  volume  difference)  and  a  energy  difference 
of  ~0.4  eV  [22].  Such  a  huge  pressure-induced  volume  collapse 
is  rather  uncommon  among  hydrides  as  well  as  inorganic  com¬ 
pounds  in  general.  However,  large  volume  collapses  (9-20%) 
under  pressure  are  observed  for  lanthanides  and  actinides  asso¬ 
ciated  with  the  valence  transition  or  localized-to-delocalized 
transition  of  f  electrons  in  these  systems.  Like  in  alkali  alumi¬ 
nohydrides,  we  observed  several  pressure-induced  structural 
transitions  for  alkaline-earth  hydrides  with  small  cation. 
For  BeH2  and  MgH2  we  have  found  four  pressure-induced 
structural  transitions  whereas  only  one  pressure-induced 
structural  transition  was  observed  for  BaH2  [49].  From  these 
systematic  studies  it  appears  that  if  the  ionic  radius  of  the  cation 
is  small,  one  can  expect  several  pressure-induced  structural 
transitions  in  main  group  metal  hydrides  and  particularly  huge 
volume  collapse  at  the  first  phase  transition.  Interestingly,  in 
alkali  aluminohydrides  we  observed  several  pressure-induced 
structural  transitions  when  they  also  have  smaller  cation 
radius. 


Fig.  6.  Total  energy  vs.  volume  curves  for  (a)  RbAlH4  and  (b)  CSAIH4  for 

different  possible  structural  arrangements.  The  arrows  indicate  the  total  energy  4.  Conclusions 
cross-over  point  where  the  structural  transitions  take  place. 

From  the  systematic  studies  on  metal  hydrides  we  have 
shown  that  one  can  pack  hydrogen  atoms  more  efficiently  in 
metal  hydrides  than  the  general  belief.  The  electron  local¬ 
ization  function  is  a  useful  tool  to  predict  site  preference 
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of  hydrogen  in  metal  hydrides.  Based  on  this  observation 
we  generalized  an  empirical  site-preference  rule  which  states 
that  “Hydrogen  prefers  to  occupy  the  interstitial  sites  where 
electrons  have  relatively  more  nonbonding  localized  nature 
than  other  possible  sites  in  metals,  alloys  and  intermetallic 
frameworks”. 

Industries  are  looking  for  hydrogen  storage  materials 
with  high  wt.%  of  hydrogen  and  good  hydrogen  absorp¬ 
tion/desorption  kinetics.  The  existing  phases  lack  one  or  both 
of  them.  One  strategy  is  to  explore  new  compounds  using 
accurate  density  functional  total  energy  calculations  and  solve 
the  crystal  structures  of  thus  identified  stable  compounds 
theoretically  and  try  to  synthesize  them  to  check  for  good 
kinetics.  We  have  made  structural  optimization  studies  on  a 
series  of  complex  hydrides  and  the  crystal  structures  of  several 
new  compounds  have  been  predicted.  As  the  known  stable 
complex  hydrides  lack  sorption  kinetics,  it  is  important  to 
identify  possible  metastable  phases  in  potential  hydrides  and 
try  to  stabilize  them  by  suitable  thermo/chemical  processing 
conditions.  From  our  high-pressure  studies  we  identified  several 
high  pressure  phases  in  these  materials  and  in  particular  we 
found  large  volume  collapse  at  the  first  phase  transition  point  in 
several  potential  hydrogen  storage  materials.  Also,  the  chemical 
bonding  behavior  of  these  high-pressure  phases  are  drastically 
different  from  that  of  ambient  pressure  phase  and  hence  it  is 
expected  that  the  metastable  phases  may  improve  the  hydrogen 
absorption/desorption  kinetics. 

Detailed  analysis  has  been  made  about  the  chemical  bond¬ 
ing  behaviors  of  complex  hydrides  using  charge  density,  charge 
transfer,  electron  localization  function,  partial  densities  of  states, 
Mulliken  population  analysis,  and  Born  effective  charges.  The 
general  feature  of  complex  hydrides  is  the  presence  of  dominant 
ionic  bonding  between  alkali/alkaline-earth  metal  and  the  B,  Al 
or  Ga  with  H  complexes.  The  bonding  interaction  between  H 
and  B,  Al  or  Ga  comprises  both  ionic  and  covalent  character. 
Generally,  the  covalent  bonding  interaction  decreases  when  one 
goes  from  B  to  Al  to  Ga  compounds. 
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